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Abstract 

We report on a search for CP asymmetry in the singly Cabibbo-suppressed decay D^ -^ 
K^K~7r~^ using a data sample of 818 pb~^ accumulated with the CLEO-c detector on the '(/;(3770) 
resonance. A Dalitz-plot analysis is used to determine the amplitudes of the intermediate states. 
We find no evidence for CP violation either in specific two-body amplitudes or integrated over the 
entire phase space. The CP asymmetry in the latter case is measured to be (— 0.03ib0.84ib0.29)%. 

PACS numbers: 13.25.Ft, 11.30.Er 



D-TCieson decays are predicted in the Standard Model (SM) to exhibit CP- violating charge 
asymmetries smaller than O{10~^) [l|]. Measurement of a CP asymmetry in the D system 
with higher rate would clearly signal new physics (NP) [2, y]. Singly Cabibbo-suppressed 
(SCS) decays via c — > uqq transitions are sensitive to NP contributions to the AC = 1 
penguin process. Interestingly, such processes do not contribute to either the Cabibbo- 
favored (c -^ sdu) or the doubly Cabibbo-suppressed (c -^ dsu) decays. Direct CP violation 
in SCS decays could arise from interference between tree and penguin processes. A non-zero 
CP asymmetry can occur if there is both a strong and weak phase difference between the 
tree and penguin processes. In charged D-meson decays, mixing effects are absent, allowing 
us to probe direct CP violation and consequently NP. 

Weak decays of D mesons are expected to be dominated by quasi two-body decays with 
resonant intermediate states. Dalitz-plot analysis techniques can be used to explore the 
resonant substructure. The intermediate structures of D^ — »■ K~^K~7t^ decay were studied 
by E687 [^ with a Dalitz-plot analysis and by FOCUS [5|] with a non-parametric technique. 
BABAR searched for direct CP asymmetries in this mode using a counting method [6|]. Using 
281 pb~^ of data, CLEO previously measured the absolute hadronic branching fractions 
and the CP asymmetries of Cabibbo-favored Z)-meson decay modes and the phase-space 
integrated asymmetry in the K'^K~7i~^ mode we study here [7|. The previous investigations 
of this decay were either limited by statistics, and did not search for CP violation, or did 
not study the resonant substructure. 

We present the results of a search for direct CP asymmetry in the decay D^ -^ K^K'n"^. 
This includes a study of the integrated decay rate, as well as decays through various inter- 
mediate states. We perform the present analysis on 818 pb~^ of e^e~ collision data collected 
at a center-of-mass energy of 3774 MeV with the CLEO-c detector |8|, y, [lOJ] at the Cor- 
nell Electron Storage Ring (CESR). The CLEO-c detector is a general purpose solenoidal 
detector that includes a tracking system for measuring momentum and specific ionization 
{dE/dx) of charged particles, a Ring Imaging Cherenkov detector (RICH) to aid in particle 
identification, and a Csl calorimeter for detection of electromagnetic showers. 

We reconstruct D~^ —>■ K~^K~7r~^, and the charge-conjugate mode D~ — > K^K~n~. 
(Charge-conjugate modes are included throughout this report unless noted otherwise.) The 
event reconstruction criteria are the same as that used in Ref. [7|]. Charged tracks are 
required to be well measured and to satisfy criteria based on the track fit quality. They 
must also be consistent with coming from the interaction point in three dimensions. Pions 
and kaons are identified using dE/dx and RICH information, when available. If either dE/dx 
or RICH information (or both) is missing we still use the track in the analysis. Detail can 
be found in Ref. [7|. We define two signal variables: 

^E = 2^ Ej — -Ebeam (1) 

i 

and 

"^BC = A'eam-IEP^P' (2) 

where Ei and p, are the energy and momentum of each D decay product, and -Ebeam is the 
energy of one of the beams. For a correct combination of particles, AE should be consistent 
with zero, and ttibc should be consistent with the D'^ mass. Fig. [1] shows AE distribution 
of data. We select candidates that have AE within ±12 MeV of zero, corresponding to 2.5 
standard deviations (a). If in any event there are multiple candidates satisfying the AE 



criterion using entirely separate combinations of tracks, we accept all of these candidates. 
Otherwise if there are multiple candidates sharing tracks we keep only the combination with 
the smallest \AE\. 
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FIG. 1: The AE distributions. Signal (|A^| < 12 MeV) and sidebands (50 MeV< \AE\ < 
100 MeV) regions are shown. 
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FIG. 2: The tjibc distributions for (a) D^ and (b) D candidates. The soUd curves show the fits 
to the data (points with error bars), while the dashed curves indicate the background. 



To determine the signal yields of the D~^ and D samples, we simultaneously fit the m-Bc 
distributions from the samples and require they have the same signal shape. For the signal, 



we use a Crystal Ball line shape function ll|, whose parameters are allowed to float. For 



the background, an ARGUS function [12[ is used with shape parameters determined from 
the events in the AE sideband (50 MeV < \AE\ < 100 MeV). We flnd 9757 ± 116 D+ 
and 9701 ± 115 D~ . Figure [2] shows the rriBC distributions of D~^ and D~ samples with fit 
functions superimposed; the total x^ is 241 for 180 degrees of freedom (d.o.f.). 

We obtain the efficiency from a GEANT-based signal Monte Carlo (MC) simulation of the 
detector. The signal MC requires one of the two D mesons in an event to decay in accordance 
with all known modes and the other one to decay to the signal mode. For the signal D meson, 
we generate events that uniformly populate phase space. The average efficiency, accounting 
for a non-uniform population density of data, is calculated as follows. The Dalitz plot of the 
data is first divided into 16 bins that are approximately equally populated. The signal yields 
are obtained from the m^c fits bin by bin and the corresponding efficiencies are calculated 
from the MC. The average efficiency is the sum of the yields divided by the sum of the 
efficiency-corrected yields. We find the efficiencies e^ for the D^ decays are (44.13 ±0.15)% 
and (43.85 ± 0.15)%, respectively. The CP asymmetry, defined as 

A - ^^/^^ - ^"/^~ C) 

"^^^ - m/e^ + N^/e-' ^^^ 

where N^ are the measured D"^ yields, is measured as 

Acp = (-0.03 ±0.84 ±0.29)% . (4) 

For the Dalitz-plot analysis, we consider the events from the signal box (jAE'l <12 MeV 
and \m^c ~ ""^d+I < 4.5 MeV/c^) corresponding to a 2.5cr range in each variable. The 
signal purity is (84.26 ± 0.10)% obtained from the mBC fit- The K~^K~7t~^ Dalitz-plot 



distribution is parameterized using the isobar model formalism described in Ref. [13|. The 
decay amplitude as a function of Dalitz-plot variables is expressed as a sum of two-body 
decay matrix elements, 

M{ml,m'i,) =^are^^''Ar{ml,m'i), (5) 

r 

where each term is parameterized with a magnitude a^ and a phase 6r for the intermediate 
resonance r, and r ranges over all resonances. We choose m'^ = m\-+^+ and m?_ = m\--^+ 
as the two independent Dalitz-plot variables. The partial amplitude Ar{rn?^,rn?_) is param- 
eterized using the Breit-Wigner sha pe with Blatt-Weisskopf form factors in the D meson 



and intermediate resonance vertices [ij], and angular dependence taken into account [13 
We use an unbinned maximum likelihood fit which maximizes the function 

^ = E 2 In £(m^„ m^_,,) - [t^] ' , (6) 

where the index i runs over all A^ events. The last term is used to constrain the signal 
fraction / to be the value /o within its error aj obtained from the m-QQ fit. The first term 
contains the likelihood function 

C(ml.ml) ^ /K.'f-)l-^l' + (1 - /) ^^»''";-'"-) . (7) 



TABLE I: Fit results for three models with different S-wave parameterizations. The A'^vr'^ S- 
wave contains contributions from 7^0(1430)*^ and a nonresonant term in fit A, from iiro(1430)'^ and 
k(800) in fit B, and from the LASS amplitude in fit C. The errors are statistical, experimental 
systematic, and decay-model systematic, respectively. 



Magnitude 



Phase 



Fit Fraction (%) 



Ko(1430)0 

ao(1450)0 
(/)(1680) 
^2 (1430)0 
NR 



Fit A [xVd.o.f. 
1 (fixed) 



= 898/708] 
(fixed) 



3.7 ±0.5" 



-0.5+1.0 



1.189 ±0.015tS&?^ 

+0.0+1.3 
0.1-0.7 
0.5+1.9 
0.4-3.6 
0.0+0.6 
0.3-0.2 



1.9 ±0.2] 

6.4 ±0.9^ 
5.1 ±0.3^ 



VQ _|_ n+6+15 

-179 ±41^15^ 

1+9 
1-15 
■0+10 
5-26 

150 ± 6lot?3 

53 ±71^11? 



25.0 ±0.6] 
12.4 ±3.3^ 



123 ± 3: 

-52 ±8 



Total Fit Fraction 



-11 

.7 ±2 



.+0.4+0.2 

-0.3-1.2 

,+3.4+7.3 

■ -0.7-5.8 

28.1±0.6t0il°i 

5.9 ± o.ituni 

0.51±0.1ll°:K?^ 

1 9 _|_ n Q+O.2+0.8 
-0.2+3.9 



14.7 ±1.^ 



-1.6-1.5 



K*' 



Fit B [x^/d.ol. 
1 (fixed) 



= 895/708] 
O(fixed) 



Ko(1430)'^ 

ao(1450)0 
(^(1680) 
^2 (1430)0 
«:(800) 



i.iOD It U.UiO_o.009-0.009 

i.5o±o.ioj:0:Ki 
i.86±o.2oJ:0:K?? 



70 ±6 



+1+16 
-6-23 



25.7±0.5l0-^l?i 
18.8±1.2t0;?l|2 

27.8 ± 0.4+0-^+°-2 



-163±3l}l^^ 
116 ± 2l|tl4 
-112 ±6+^+}9 0.51 ±0.1llO-Oi+0-?I 



4.6 ±0.6" 



0.3-0.4 
0.5+7.2 
0.3-1.. 



7.6 ±0.8' 



-0.5+2.4 
-0.6-4.8 



171 ± 4; 

2.30±0.13to:?i-o:29 -87 ±6 

Total Fit Fraction = (86.1 ±1.1 



4-12 
0+24 
2-11 
-2+15 
-3-10 



1.7 ±0.4 
7.0 ±0.8' 



-0.04-0.15 
+0.3+1.2 
-0.2-0.7 
-0.0+3.5 
-0.6-1.9 



LASS 

oo(1450)0 
(?!)(1680) 
K; (1430)0 



Fit C [xVd.o.f. 
1 (fixed) 
3.81 ±0.06: 
1.193 ±0.015" 



,+0.05+0.13 
-0.05-0. 46 

-+0. 003+0.021 
-0.010-0.011 

+0.14+0.68 
-0.03-0.38 



912/710] 
O(fixed) 25.3 ± 0.5 

25.1±2+^t^ 40.6 ± o.8+0:^j:9:? 



+0.2+0.2 
-0.4-0.7 



176 ± 21^1^ 



28.6±0.4+0:210;2 



1.73 ±0.07: 

1 71 +n lfi+O-02+0.41 
i./i It U.iD„o.02-0.77 



4.9 ±0.7; 



-0.1+2.2 
-0.4-2.3 



Total Fit Fraction 



122 ± 2; 

-72 ±8 
146 ± 9; 
= (101.5 ±0 



1+8 

1-10 

-2+10 

-2-22 

0+34 

7-11 



6.0 ±0.4- 



-0.9+5.5 
-0.2-2.4 



0.42±0.08tO:Kl^ 



0.7 ±0.2: 



-0.0+0.7 
-0.1-0.3 



where 



and 



A/Lig = e{m'l,m'i)\M\'^dm'ldm'^_ 
Abg = / F}^g{rri^,rri^_) d7n'^d7n'i_ 



"+' 



(8) 



(9) 



are the normalization factors, and e{m^,m'^) and Fbg(m^,?72^) are efficiency and back- 
ground functions. The fit parameters are a^, 0r and /. 

We determine the efficiency £:(m^, m^) using the same signal MC sample described before. 
The efficiency function is parameterized by a cubic polynomial in (m^, m^) multiphed by 
threshold factors T{'m'^_^^^^ - m\]pr,y) x T{m'i^^^ - mi^p^y) x T{zmax - z;pz), where 



T{x;p) 



sm{px) , < px < 7r/2 
1 , otherwise 



(10) 



z = m|.+j^_, m^^^^ or Zmax is the maximum value of m| or z in this decay, "Pxy and p^ 
are the fit parameters. The threshold factors are used to account for tracking inefficiency 
at the Dalitz-plot corners, where one of three particles might be produced with very low 
momentum and escape detection. 

Figure [2] shows that the background is significant. To construct a model of the back- 
ground shape Fbg(m^,m?_), we select events from the sideband region (24 < \I^E\ < 42 
MeV and |mBC ~ '^^d+\ < 9 MeV/c^). There are 12324 events, about 3.5 times the amount 
of background we estimate in the signal region, which is dominated by random combinations 
of unrelated tracks. Although the background includes (p and K* mesons combined with 
random tracks, these events will not interfere with each other. Thus the shape is param- 
eterized by a two-dimensional quadratic polynomial with terms representing non-coherent 
contributions from and K* meson decays, multiplied by the threshold factors. 
_We consider_fifteen intermediate states, 07r+, 0(168O)7r+, K*^K+, F*(1430)0K+, 
K*{U10fK+, K;{im)'K+, k{800)K+, /o(980)7r+, /o(1370)7r+, /o(1500)7r+, /2(1270)7r+, 
/^(1525)7r+, ao(980)°7r+, ao(1450)°7r+ and a2(1320)°7r+, as well as a nonresonant (NR) con- 
tribution. The parameters of the established resonances are taken from Ref. [l5l| . except for 
the /o(980) which is taken from Ref. [l6| and the ao(980) taken from Ref. [17|. A complex 
pol e function is used to model the k(800) with pole position at s^ = (0.71 — iO.Siy GeV^ 
[18|. The nonresonant contribution is modeled as a uniform distribution over the allowed 
phase space. For the K~7r~^ S'-wave states in the decays, we also consider the LASS am- 
plitude as described in Ref. 19|, l20|, instead of a coherent sum of the states Kq{14:30)^K~^, 



k{800)K~^ and the nonresonant term. 

This study is sensitive only to relative phases and magnitudes. The mode K K'^ is 
assigned to have zero phase and unit magnitude. We choose the same phase conventions for 
the intermediate resonances as E687 [J] used. 

We begin to fit the data by considering only the three components K , 0, and J'i'Q(1430)° 
and obtain a result consistent with E687. To present a relative goodness-of-fit estimator, 
we divide the Dalitz-plot region into bins with dimensions 0.05 (GeV/c^)^ x 0.05 (GeV/c^)^ 
and calculate x^ as 

721 

X^ = -2E-an(^), (11) 



where Ui (pi) is the observed (expected) number of events in the ith bin [15|]. We find 
X^ = 1292 for (721 — 5) d.o.f. in the "three resonances" fit, where 721 is the number of valid 
bins inside the kinematically allowed region. 

Our twenty times larger statistics than E687 require a better model. We determine 
which additional resonances to include by the following procedure: starting from the three 
resonances and adding new resonances one at a time, we choose the best additional one at 
each iteration, stopping when no additional resonances have fit fractions (FF) more than 3a 
from zero. The fit fraction is defined as 

|ar^r|^ dm\dm^_ 
FFr = ^, . (12) 



The results of our fits are presented in Table [H We find that three fits (denoted as A-C) 
describe the data with similar quality. The only difference among them is in description 
of the K~Ti^ 5'-wave contribution, which is represented by the ^^0(1430)° and NR in fit A, 
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FIG. 3: (a) The Dalitz plot for D'^ -^ K^K vr"*^ candidates, (b)-(d) Projections of the results of 
the fit B (line) and the data (points). The dashed line shows the background contribution. 



by A'olWSO)'^ and k(800) in fit B, and by the LASS amplitude in fit C. Figure [3] shows the 
Dalitz plot for the D^ -^ K^K~'k^ candidates and three projections of the data with the 
result of fit B superimposed. 

We generate seven sets of GEANT-based signal MC samples with the model from fit A. 
Each set contains about the same size as in the data. We find that the fits can recover the 
input magnitudes and phases within their errors. 

Fit B gives the best agreement with the data; thus we choose it to search for CP- violation 
{CPV). The resonances in D^ {D~) decays are allowed to have different magnitudes, ar + K 
{ttr — K), and phases, 6r + 0r (^r — (pr), in the decay amplitude M. {^A). We perform a 
simultaneous fit to D^ and D~ samples. In the fit, the signal term in Eq. ([7]) is replaced by 



C 



fe+{mi,m'i)\M\' 



Slg 



e'^ (m^, m'i)\M\'^ dm^dml 



(13) 



for the D^ sample and by 

T. = 



£,,, = fe~i<,rnl)\M\' ^^^^ 



for the D sample, where e are efficiency functions obtained from the D signal MC 

separately. We cannot determine the relative magnitude and phase between D+ and D~ 

— *o 
directly, and assume 6 = and </> = for the K resonance. The free parameters in the fit 

are br/ar, ar, Sr, (t>r and /. 



Following Ref. J21|, we also compute the CP-conserving fit fraction as 



/ \2arAr\'^ dm'^ dw? 
j{\M\'^+\M\^)dmldm^. 



mcpc)r = ,.,.:,/: ,;^,..L. - 2 > (15) 



the CPV fit fraction as 



^ ^' J{\M\^ + \M\')dmldmV ^ ' 

and the CPV interference fraction (IF) as 

/ ^pafcC*''''- cos(0fc - (j)r)Ak] hrAl dm\ dml 



IF 



ky^r 



r 



J{\M\'^ + \M\'^) dmldm'^_ 



(17) 



The CP-conserving fit fraction is the same for the D^ and D~ by construction. The CPV 
fit fraction defined by Eq. ( IT6l) is sensitive to CP violation in the resonant decay. The CPV 
interference fractions of Eq. flT7|) sum over the contribution proportional to a^e'^'^^'^hr so they 
are sensitive to CP violation in interference between resonances. The phases are important 
and allow the possibility of cancelation in this sum. 

In Table [TTl we report the magnitude asymmetries hr/ar-, phase differences ipr and fit 
fraction asymmetries. The fit fraction asymmetry is computed as the difference between 
the D^ and D~ fit fractions divided by the sum. The largest fit fraction asymmetry, for 
the A'2(1430)'^, is l.To", and occurs because the fit fraction for the -R'2(1430)° is small. The 
CP-conserving fit fractions and the 95% confidence level (CL.) upper limits for CPV fit 
fraction, CPV interference fraction, and the ratio of CPV interference to CP-conserving fit 
fraction are given in Table IIIII We notice that the CP-conserving fit fractions are consistent 
with those of fit B in Tabled Figure H] shows the difference of the Dalitz-plot projections of 
data and fit between D^ and D~ decays. 

We calculate an integrated CP asymmetry across the Dalitz plot, defined as 

•^"'' = / iSrJw *"- *"-/ / *"+ *"-■ (18) 

We obtain Acp = (—0.4 ± 2.01o;51q;3)%, where the errors are statistical, experimental sys- 
tematic, and decay-model systematic, respectively. 

Using the same counting technique as in Ref. |6|, we examine CP asymmetries {Acp) in 
the and K regions by requiring the K^K~ and K~t:^ invariant mass to be within 15 



TABLE II: The magnitude asymmetries h^/aj., phase differences (p^ and asymmetries on the D'^ and 
D^ fit fractions from fit B. The errors are statistical, experimental systematic, and decay-model 
systematic, respectively. 



r 


h/a (%) 


4>{°) 


FF asymmetry(%) 


^*u 


(fixed) 


(fixed) 


-o.4±2.ot°:^J:°J 


Kliimr 


A±2,tltl 


-l±6+o+? 


8 ± Qt\t\ 


</. 


-^.7±l2tlltr2 


3 ± 3+?+? 


1 Q _|_ 1 f:+0.0+0.2 
1.8 It i.D_Q4„o.l 


ao(1450)0 


-10±7±2t^ 


4 ± ^tltl 


-19 ±12+^+^1 


(/>(1680) 


-4±llt^+6 


3 ± 6 ± 2l^ 


-9 ± 221^^+9 


i^2 (1430)0 


90+12+1+3 

^"3-11-7-7 


c+5+1+3 

-^-4-3-1 


43±19l}+_'i2 


k(800) 


-6±6l?t^ 


^±Qttt\ 


-12 ±11^0^1^ 



TABLE III: The CP-conserving fit fractions from Eq. (jlSp and the 95% confidence level (C.L.) 
upper limits for CPV fit fraction from Eq. (J16p . CPV interference fraction from Eq. (J17p . and 
the ratio of CPV interference to CP-conserving fit fraction. The 95% C.L. upper limits include 
statistical and systematic effects. 



FF(CPy) IF Ratio 
(xlO-3) (xlO-3) (%) 
Component FF(CPC)(%) (95% C.L. upper limits) 



K*' 


25.7 ±0.5 


(fixed) 


O(fixed) O(fixed) 


KliUSOf 


18.8 ±1.2 


<4.3 


< 12.6 


<8.5 


<P 


27.8 ±0.4 


<0.6 


<0.5 


<0.17 


ao(1450)° 


4.7 ±0.6 


< 10.8 


<31.6 


<45 


(/)(1680) 


0.50 ±0.11 


<0.9 


<4.6 


<89 


Kl{U30)^ 


1.8 ±0.4 


<6.9 


<3.9 


<22 


At(800) 


7.0 ±0.8 


<4.2 


< 17.2 


<25 



and 10 MeV/c^ of the nominal and K*° masses ^. We find Acp (-0.9 ± 1.4 ± 0.7)% 
and (0.3 ± 1.8 ± 0.6)% for the (p and K region, respectively. 

Systematic uncertainties from experimental sources and from the decay model are con- 
sidered separately. Our general procedure is to change some aspect of our fit and interpret 
the change in the values of the magnitudes, phases, fit fractions, br/ar, (pr, and fit fraction 
asymmetries as an estimation of the systematic uncertainty. 

Contributions to the experimental systematic uncertainties arise from our model of the 
background, the efficiency and the event selection. Our nominal fit fixes the coefficients of 
the background determined from a sideband region. To estimate the systematic uncertainty 
on this background shape, a fit is done with the coefficients allowed to float and constrained 
by the covariance matrix obtained from the background fit. Similarly, to estimate the sys- 
tematic uncertainty on the efficiency parameters, we perform a ffi with the coefficients of 
efficiency allowed to float constrained by their covariance matrix. To estimate the systematic 
uncertainty on MC simulation for the particle identiflcation, a flt is done with new efficiency 
parameters obtained from the weighted MC sample by the efficiency ratios of data to MC 
depending on each particle's momentum. To estimate the event selection uncertainty, we 
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FIG. 4: The difference of the Dahtz-plot projections of data (points) and fit (fine) between D"*" 
and D" decays. 



change the AE and m-Qc selection criteria in the analysis. These variations to the standard 
fit are the largest contribution to our experimental systematic errors. In the CP asymmetry 
search, we take the background fractions and shapes to be the same for the D~^ and D~ sam- 
ples. To estimate the uncertainty on the supposition, we perform a fit with the background 
determined separately. 

The systematic error due to our choice of D~^ -^ K~^K~7r~^ decay model is evaluated as 
follows. We change the standard values of the radial parameter in the Blatt-Weisskopf form 
factors [ij] for the intermediate resonance decay vertex (1.5 GeV~^) and the D~^ vertex 
(5 GeV~^) both to 1 GeV~^. Fits with constant width in the Breit-Wigner functions are 
considered. To compute the uncertainty arising from our choice of resonances included in 
the fit, we compare the result of our standard fit to a series of fits where each of the reso- 
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nances, K*{U10f, /o(980), /o(1370), /o(1500), /2(1270), /2(1525), ao(980)° and a2(1320)°, 
is included one at a time. These variations to the standard fit result in the largest contri- 
bution to systematic errors associated with our decay model. The masses and widths of the 



intermediate resonances are allowed to vary within their known uncertainties [15|. For fit 
C, we vary the parameters in the LASS amplitude within their uncertainties. 

We take the maximum variation of the magnitudes, phases, and fit fractions, br/cir, 4>r, 
and fit fraction asymmetries from the nominal fit compared to the results in this series of 
fits as a measure of the experimental systematic and decay-model systematic uncertainty. 
Table HVl shows the systematic checks on the integrated CP asymmetry defined in Eq. (fTSll . 
Apart from the sources discussed above, we also consider different models from fit A or C; 
the variations are small. 

TABLE IV: Sources contributing to systematic uncertainties on the integrated CP asymmetry 
defined in Eq. ([IS]). 



Source 


Variation (%) 


Background shape 
Efficiency parameters 
Particle identification 


-0.01 
+0.02 
+0.06 


Event selection criteria 


+0.18 


Background (in) dependent fit 


-0.52 


Form factors 


+0.21 


Width parameterization 
Choice of resonances 
Resonant masses and widths 


-0.15 

+0.61 
-0.33 
+0.09 
-0.08 


Fit A 


+0.07 


Fit C 


-0.15 



We estimate the systematic uncertainty on the CP asymmetry defined in Eq. (jlj). The 
contributions from various identified sources are listed in Table |Vl The uncertainty due to 
selection criteria is estimated by doubling the AE signal window. We evaluate an uncertainty 
for the background shape by floating its parameters in the fit instead of fixing them from the 
values obtained form the AE sideband. We use the CP-conserved channels D~^ — > K^ti^ti^ 
and D^ -^ K'tt^tt^ as control modes to assign the systematic uncertainty on MC simulation 
due to possible efficiency difference on positive and negative charged kaons and pious. 

TABLE V: Systematic uncertainties on the CP asymmetry defined in Eq. ^. 



Source Variation (%) 

Selection criteria ±0.25 

Background shape ±0.02 

MC simulation ±0.15 

Total ±0.29 



In conclusion, we have analyzed the resonant substructure in D'^ -^ K^K vr^ decay and 
searched for CP violation in the decay and its intermediate resonances. We measure the 
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overall CP asymmetry in D^ -^ K+R-tt^ decays to be (-0.03 ± 0.84 ± 0.29)%. The limit 
is more restrictive than the one found previously by BABAR J6|. We use five resonances and 
K~7r~^ S'-wave states to model the Dalitz plot with results shown in Table HI The K~7r~^ 
5- wave can be equally well described by a coherent sum of Kq{14:30)^ and nonresonant 
amplitude or ^0(1430)° and /t:(800), or the LASS amplitude. Choosing the second model 
we measure the CP asymmetries for all submodes, shown in Table HTl and lllli The measured 
CP asymmetries are consistent with the absence of CP violation. We find Acp defined 
in Eq. flTSl) to be (—0.4 ± 2.0to;5lo!3)%- The Acp is sensitive to an asymmetry in shape 
between the D~^ and D~ samples, but does not depend on their yields. 
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